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ABSTRACT 


The  temperature  dependence  of  electron  mobility  of 
high  purity  (n  < (T~x  fo^  cmi”^)  LPE  Ga^_J>Al^s  layers 
(x  < .18)  has  been  studied  in  the  range  25 °K  - 110°K  to 
measure  the  extents  of  alloy-like  and  ionized  impurity 
scattering.  Values  of  the  alloy  scattering  parameter  /e 
(which  may  include  a contribution  from  space  charge  scat 
tering)  determined  for  these  samples  were  found  to  lie 
in  the  range  .36  eV  to  .51  eV,  the  average  being  .44  eV. 
The  compensation  ratio  was  found  to  be  about  2,  and 
independent  of  x.  D 
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I . INTRODUCTION 
A.  Alloy  Scattering 


Alloy  scattering  of  electrons  in  a semiconductor 
alloy  such  as  Ga^  xAlxAs  is  caused  by  randomness  in  the 
distribution  of  the  alloying  constituents.  An  expression 


for  the  alloy  scattering  component  of  mobility  was 

Cl  2} 

first  derived  by  Brooks v ’ as 


(2tt) 


q 'fi4  N, 


* (kT)^  m*5"2  x(l-x)  (AE)  2 


(1) 


where  q is  the  electronic  charge,  is  (Plank's  constant)  / 2tt, 

N is  the  density  of  alloying  sites , k is  Boltzmann  s 
d 

constant,  T is  the  absolute  temperature,  m*  is  the  effec- 
tive electron  mass,  x and  (1-x)  are  the  mole  fractions 
of  the  binary  end  compounds,  and  AE  is  an  energy  step 
which  was  assumed  equal  to  the  band-gap  difference  between 
the  end  compounds . 

f 3) 

In  a recent  paper,  ' Harrison  and  Hauser  derived  a 
similar  expression  for  alloy  scattering.  They  assumed  the 
crystal  to  have  a uniform  background  potential  with  square 
wells  of  depth  AE  at  random  sites  associated  with  one  of 
the  alloying  constituents.  The  square  well  potential  was 
assumed  to  extend  over  a spherical  region  of  radius  equal 
to  the  nearest  neighbour  separation,  a choice  they  ad- 
mitted as  being  somewhat  arbitrary.  They  obtained  the 


2) 


energy  dependent  scattering  time  xA(e)  as 


ta<£> 


N e ^ 
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3/2  tt  (kT)*  m*  2 x(l-x)  (AE) 


where  e = (E-Ec) /kT  is  the  electronic  kinetic  energy 

normalized  to  kT.  Multiplying  Eq.  (2)  by  — %■  •a—  gives 

3/tt  m* 

an  expression  for  vA  that  is  identical  to  Brooks'  expres- 
sion except  for  a constant  factor.  Thus,  for  a given  set 
of  physical  parameters,  (uA)gr^  = 1.85  (^A)HH,  where  Brk 
and  HH  denote  the  use  of  expressions  (1)  and  (2)  respec- 
tively to  calculate  uA- 

In  Brooks'  original  formulation,  AE  was  taken  to 

(4) 

be  the  band  gap  difference  between  the  two  end  compounds.  ' 
For  ternary  alloys  where  both  end  compounds  have  direct 
gaps,  (such  as  Ga^_xInxAs) , there  is  no  ambiguity.  How- 
ever, if  one  of  the  end  compounds  is  direct  and  the  other 
is  indirect,  as  in  Ga, _ A1  As,  then  the  question  arises 
whether  AE  is  the  band  gap  difference  (AE  =0.73  eV)  or 

the  direct  band  gap  difference  (A(Epc-Erv)  = 1.52  eV) . 

(3) 

Harrison  and  Hauserv  J have  suggested  that  instead  of  AE  , 

6 

AE  should  be  taken  as  the  conduction  band  edge  discon- 
tinuity AEc  between  the  end  compounds.  Again,  taking  AE 
to  be  A(EC  - Eyac)  gives  a value  0.55  eV,  while  A(Ep  c - Eva<;) 
gives  1.34  eV.  Finally,  Ferry ^ has  used  the 


3) 


electronegativity  differences  to  predict  that  AE  for 
Ga,  Al  As  should  be  , 12  eV  at  x = .5. 

1 " a X 

It  is  clear  that  depending  on  which  of  the  several 
"justifiable"  values  of  AE  (ranging  from  . 12  eV  to  1.52  eV) 
is  selected,  and  which  of  the  two  expressions  (1)  and  (2) 
is  used,  alloy  scattering  in  Ga^ 
as  being  negligible,  moderate,  or  dominant.  Thus  until  a 
clearer  theory  of  alloy  scattering  is  developed,  we  must 
turn  to  experiment  to  determine  the  value  of  AE. 

B.  Space  Charge  Scattering 

Space  charge  scattering  in  semiconductors  is  caused 
by  space  charge  regions  that  form  around  localized  compo- 
sition and  impurity  concentration  inhomogeneities.  This 
mechanism  was  first  discussed  by  Weisberg^  to  explain 
anamolously  low  mobilities  in  GaAs . Later,  Conwell  and 
Vassel^^  proposed  a simple  model  where  the  space  charge 
scattering  time  was  given  by 

tsc(e)  = (NsQv)'1  (3) 

from  which  the  space  charge  mobility  component  is  obtained 
as 

ugc  = 2.4xl09/(m*/mo)%  T%  NSQ  cm2/Vsec  (4) 

_ 3 

where  Ns  is  the  volume  density  (cm  ) of  the  space  charge 
scattering  centers  and  can  vary  widely  depending  on  the 


^lxAs  can  be  represented 


4) 


quality  of  the  sample,  Q is  the  cross  sectional  area  of 
2 

the  centers  (cm  ) and  v is  the  electron  velocity  (cm/ sec) . 

Subsequently  Katoda  et  al  and  others (9.10)  ^ave 
used  Eqs.  (3)  and  (4),  which  have  the  same  temperature 
and  energy  dependence  as  alloy  scattering,  to  describe  space 
charge  scattering.  However  Weisberg^* 

dependence  of  the  scattering  cross  section  Q to  be  (n/T) 


obtained  the  temperature 

1/3 


so  that 


usc  * n 


(5) 


n being  the  free  carrier  density. 

We  are  inclined  a priori  to  agree  with  Eq.  (5). 

However,  in  this  work  we  found  it  necessary  to  ignore  the 

(n/T)  3 dependence  of  Q.  Since  we  were  attempting  to 

measure  the  extent  of  alloy  scattering  experimentally, 

merging  space  charge  scattering  with  alloy  scattering 

would  simplify  our  analysis  considerably.  Furthermore, 

our  results  showed  that  space  charge  scattering  does  not 

dominate  over  alloy  scattering  in  our  high  pruity  samples, 

(9) 

as  it  does  in  samples  studied  by  Kaneko  et  al.  ' and 
Stringfellow.  This  offered  a post-priori  justification 

for  making  this  approximation. 

We  use  the  term  "alloy-like"  scattering  to  collectively 

_ v 

denote  the  scattering  mechanisms  with  the  (eT)  ^ relaxa- 
tion time  dependence.  We  define  the  "alloy-like"  scatter- 
ing parameter,  Eg,  of  a given  ternary  sample  to  be  the 


value  of  AE  in  Eq.  (2)  required  to  fit  the  experimental 
y vs  T data,  when  ionized  impurity  and  lattice  scattering 
is  also  considered.  Thus  defined,  Eg  absorbs  any  arbi- 
trariness in  the  constant  factor  in  Eq.  (2),  and  includes 
contributions  (if  any)  from  space  charge  scattering.  Our 
choice  of  Eq.  (2)  over  Eq.  (1)  in  defining  Eg  rests  entirely 
on  the  fact  that  the  derivation  of  the  latter  is  unpublished. 
Using  Eq.  (1)  instead  of  Eq.  (2)  to  define  Eg  would  give 
values  higher  by  a factor  1.36. 

In  this  study  we  report  on  the  experimental  deter- 
mination of  Eg  for  nine  samples  of  high  purity  Ga^_xAlxAs 
(x  < .18)  grown  by  LPE. 

II . EXPERIMENT 

A.  Principle  of  the  Experiment 

The  alloy  scattering  parameter  Eg  can  be  determined 
experimentally  by  separating  the  various  scattering  mech- 
anisms in  Ga^_xAlxAs.  The  Hall  mobility  y^  of  a ternary 
alloy  is  a function  of  x,  T,  B,  and  Eg,  where  B is 
the  magnetic  field  and  Ng  * (N-^  + NA)  times  the  screening 
factor.  By  estimating  the  various  physical  parameters, 
for  a given  x and  T (as  described  in  section  III-B)  the 
polar  optical,  deformation  potential  and  piezoelectric 
scattering  time  constants  can  be  determined.  The  ionized 
impurity  scattering  relaxation  time  can  be  determined 


as  a function  of  N , while  alloy  scattering  can  be  ex- 

S 

pressed  in  terms  of  Eg,  both  of  which  are  unknown.  Thus 
the  net  relaxation  time,  the  drift  mobility,  and  the  Hall 
mobility  of  Ga,  A1  As  can  be  calculated  as  functions  of 
Ng  and  Eg,  as  described  in  section  III-C. 

To  measure  Ng  and  Eg,  use  is  made  of  the  fact  that 
varies  as  T 2 while  xA  varies  as  T . By  measuring 
over  a temperature  range,  the  alloy  and  the  ionized 
impurity  components  can  be  separated  from  the  residual 
scattering  (i.e.,  after  phonon  scattering  has  been  sub- 
tracted from  the  total) . 

B.  Experiment 

The  Hall  mobility  measurements  were  made  on  ten 
layers  of  high  purity  n Ga,  A1  As  (x  < .18)  (and  one 
reference  layer  of  undoped  n GaAs) , grown  by  LPE  at  700  - 
680°C.^^  The  layers  were  9-14  microns  thick  and  were 
grown  on  nonconverting  S.I.  GaAs  substrates.  Two  layers 
were  grown  from  each  of  five  Ga-Al-As  melts  of  increasing 
A1  concentration.  The  first  layers  will  be  referred  to 
at  type-A,  and  the  second  as  type-B  layers.  Van  der  Pauw 
Hall  samples  were  fabricated  from  each  layer,  and  tin  dot 
contacts  were  alloyed  under  hydrogen  using  a strip  heater. 
Good  ohmic  contacts  were  obtained.  Van  der  Pauw  measure- 
ments were  made  on  samples  at  B = 2 KGauss,  at  5°  intervals 
from  25°K  to  110°K. 


The  Hall  mobility  and  the  Hall  carrier  density 
per  unit  area  QH  obtained  from  the  measurements  are  plotted 
against  temperature  in  Figs.  1 and  2 respectively.  Room 
temperature  mobilities  were  also  measured  and  are  reported 
elsewhere.  The  layer  compositions  were  estimated  from 

5°K  photoluminescence  measurements.  Although  a SIMS  analy- 
sis on  a test  layer  indicated  a 127„  increase  in  x in  going 
from  the  surface  to  the  substrate  interface,  we  assumed  for 
simplicity  that  x was  constant  and  equal  to  the  value  mea- 
sured by  p.l.  The  layer  thicknesses  were  measured  by  cleave 
and  stain,  and  are  listed  along  with  the  p.l.  data  in  Table 
A.  The  last  column  lists  x as  estimated  from  theoretical 
solidus  data  for  the  Ga-Al-As  system^^  ’ , and  shows 

good  agreement  with  the  experimental  values . 

III.  ANALYSIS 
A.  Analysis  Technique 

The  data  was  analyzed  by  computer.  For  each  pair  of 
values  of  x and  T,  the  following  calculations  were  performed. 
(1)  The  various  physical  parameters,  namely  the  effective 
electronic  mass  (m*) , the  static  and  optical  dielectric 
constants  (K  , K ),  the  longitudinal  and  transverse  elastic 

o u 

constants  (c^,  ct),  the  deformation  potential  (Ec)  , the 
piezoelectric  constant  (h^)  and  the  longitudinal  optical 
phonon  energy  (LOPE,  hw^),  were  estimated,  and  used  to 
calculate  the  three  lattice  scattering  times  as  functions 


(2)  A subroutine  "MOBILITY"  was  used  for  calculating  y^ 

2 

and  y^,  given  x,  T,  B,  EBSQ  (test  value  of  Eg)  and  NS 
(test  value  of  Ng) . 

2 

(3)  For  a series  of  test  values  of  Eg,  the  parameter  NS 
was  adjusted  to  generate  a y^  equal  to  the  measured  value. 
The  corresponding  drift  mobility  y^  and  Hall  factor 

rK  = were  a^so  calculated. 

(13) 

(4)  The  free  carrier  density  was  calculated  from' 

„ , rHQH 
sc 

where  d is  the  layer  thickness  and  d is  the  sum  thickness 

b C 

(13) 

of  the  surface  and  interface  depletion  regions.  ' The 

surface  and  interface  barrier  heights  used  for  estimating 

(13") 

dgc  were,  respectively'  ' 

SBH  - (0.6  + 0 . 8x)  eV  (7) 

IBH  = (0.75  + l.lx)  eV  (8) 

The  variation  of  SBH  with  x was  estimated  from  reference 

(14) , while  for  IBH  it  was  assumed  equal  to  the  variation 

of  AE  with  x. 
c 

(5)  From  the  values  of  NS  and  N calculated  in  the  previous 
two  steps,  ND  + NA  and  hence  ND  and  NA  were  calculated 
using  Eqs.  (24)  and  (25). 

(6)  Thus  for  each  test  value  of  EBSQ,  values  of  NA  were 
generated  at  each  temperature.  Since  in  n-type  material, 

"NS,  ND,  NA,  N and  EBSQ  represent  test  values  of  the 
physical  parameters  Ng,  NJ,  N^,  n and  Eg. 


9) 


* 
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the  ionized  acceptor  density  is  temperature  independent, 


the  constancy  of  NA  with  T was  used  as  the  criterion  for 

selecting  Eg  from  the  series  of  test  values  of  EBSQ.  The 

corresponding  values  of  NA,  ND  and  N were  taken  to  represent 
★ 

Na,  Nd  and  n. 

(7)  To  summarize,  the  program  took  as  its  input  the  values 

2 

of  x,  d,  B,  sets  of  (T,  yR,  QR) , and  test  values  of  Efi , and 

printed  out  values  of  NA,  ND  arid  N vs  T for  each  test 
2 

value  of  Eg. 

B,  Estimation  of  Parameters 

The  relative  effective  mass  (m*/mQ)  for  electrons 

was  assumed  to  be  .067  in  GaAs^^  and  .014  in  A£.As^  ^ 

at  temperatures  below  110°K.  Its  variation  with  x was 
/ 1 8 \ 

calculated  from'  ' 


(?) 
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X 


+ i^T> 


-1 


(9) 


The  static  and  optical  dielectric  constants  for 

Ga,  Ail  As,  denoted  by  Kx  and  Kx  respectively,  were 
i-x  X so 

calculated  from  the  expression 


(18) 


Kx-1 

Kx+2 


,K^-1 
K°+2 


K1-! 

■)  (1-x)  + (V^)x 
Ki+2 


(10) 


where  K°  and  K1  refer  to  values  for  x = 0 and  1 respectively. 

K°  was  obtained  as^^ 

3 , 

K°  = 12.91  (1  - (300-T)  x 10_H}  , (11) 

s 

while  K° , K1,  and  K1  were  assumed  to  remain  proportional 
o s o 

to  K°  at  all  temperatures,  their  room  temperature  values 

(19  20} 

being  10.91,  10.9  and  8.5  respectively v ’ . 


*See  footnote  on  page  8 
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(21) 

According  to  Keyes'  ' , the  elastic  constants  of 

III-V  compounds  /ary  only  with  the  lattice  parameter  a 

which  for  Ga,  AJ l As  is  virtually  independent  of  x.  The 

(22) 

longitudinal  and  transverse  elastic  constants  c^  and  ct 
of  GaAs  were  therefore  used  for  (Ga,A£)As,  so  that  ^0) 

cz  = 1.444  x 1012  (1-5 . 4x10_5T)  dynes/cm2  (12) 
and  = .4905  x 10^2  dynes/cm2  (13) 


(23) 


, dEc/dP  (where  P is  the 


Also  according  to  Keyes 
hydrostatic  pressure)  is  the  same  for  all  III-V  compounds 
with  the  zinc  blend  structure.  Since  the  elastic  constants 
are  assumed  to  be  independent  of  x,  it  follows  that  the 
variation  of  Ec  with  strain  is  also  x- independent . Tie 
deformation  potential  of  the  r band  edge  in  Ga,  Ai.  As 

C i “ X X 

was  therefore  assumed  constant  at  7.0  eV,  the  value 

(22) 

determined  by  Wolfe  et  al'  ' for  GaAs. 

Estimation  of  the  piezoelectric  constant  h-^^(V/cm) 
was  based  on  the  assumptions  that  (i)  the  dimensionless 
product  es^x4//e14  ^es  = stat^c  permitivity)  Is  constant  for  all 
III-V  semiconductors,  and  (ii)  e^,  which  is  the  piezoelectric  charge  per 
unit  strain,  varies  linearly  with  x.  The  first  assumption  was  verified 
by  testing  measured  values  of  h^  and  e^  for  five  direct 
band  III-V  compounds ^2^  , for  which  e h^/e^  was  found 
to  be  constant  within  570.  h-^^(x)  was  obtained  from 


h14  H h14(x) 


xKah14  + (1-x)  KshU 


K. 


(14) 


11) 


where  h°^  was  assumed  to  be  1.4  x 10  V/ cm  , and  h-^ 
was  obtained  as  =1  x V/cm  by  assuming  e-^(A2As)  to 

C25) 

be  .1,  the  value  for  e-^(GaP) , as  estimated  by  Hauser' 

The  temperature  variation  of  the  longitudinal  optical 

phonon  energy  (LOPE,  hw.  ) for  GaAs  was  obtained  from  the 

o 

consideration  that 


- <VM) 


1/2 


where  M is  the  reduced  mass , and  that 
dc , 


1_ 

CJc 

so  that 


3T 


'l  ~ 


5.4  x 10 


-5 


dcj 


2,0 


2o 


dT 


s - 2.7  x 10 


-5 


(15) 


The  variation  of  hu>0  with  x for  the  GaAs -like  and  the 
A2As-like  phonon  branches,  denoted  by  L0PE(1)  and  L0PE(2) 
respectively,  in  the  range  x < .3,  was  obtained  from  the 
data  of  Ilegems^0>26)  . The  ^ase  temperature  values  were 
obtained  from  references  (22)  and  (26)  respectively,  to 
give  the  expressions 

L0PE(1)  = 36.465  [1  + 2.7xl0-5  (77-T)]  - .625x  meV 


and 


(16) 


LOPE (2)  = 44.7  [1  + 2.7xl0'5(300-T)]  + 15. 5x2  meV 


(17) 


C.  Calculation  of  Relaxation  Times  and  Mobilities 
The  expressions  for  the  deformation  potential  and 
the  piezoelectric  scattering  relaxation  times  were  taken 


12) 


from  reference  (22): 

TDp(£)"1  = 4 . 167xl019  E2  c’1  (m*/mo)-3/2  T3/2e1/2sec~1 

(18) 

f N-l  1 in7u2  r4  .3  vm*,l/2Tl/2  -1/2  -1 

TpE(e)  = 1.0524x10  h14L  — + - J (— ) T e sec 

C 1 (19) 

Polar  optical  scattering  for  the  binary  compounds  was 
represented  by  a time  constant v ' 


'rpo(e) 


-1 


T?0  + 


(20) 


'e-z 


i 


l 


where  z ■ 
e>z  only. 


(hui^/kT)  , and  the  second  term  is  considered  for 


-po 


was  calculated  by  obtaining  an  expression 


for  UpQ  in  terms  of  XpQ  based  on  Eq.  (20),  and  equating 

(27) 

it  to  the  value  of  vipn  obtained  fromv 

ru  i i 

(—  - — ) 
vx  v* 

, _2 ,,,  \ .2357  o ^ 

“P0<C"  /VSeC>  ' [h„„(eV)]^  ^T?71  ^ 


3 ^ (ez-l) 


(21) 

where  x(z)  can  be  approximated  by  3/ttz/8  for  z > 5^22^  , 
a condition  that  is  valid  for  GaAs  below  85°K. 

The  presence  of  two  optical  phonon  branches  in 
Ga, _ A Z As  was  treated  by  assuming  that  the  densities  of 
GaAs-like  and  A£As-like  phonons  are  proportional  to  (1-x) 
and  x respectively.  Thus  polar  optical  scattering  in 
Ga^  xA£xAs  was  represented  as 


Tp0(e) 


•1  = 1-x  1_ 

/e 

TP0 


'e-z. 


x 

rP0 


l-  + £- 

/c 


/e-z. 


(22) 


13) 


where  z±  = LOPE(l)/kT,  z2  - L0PE(2)/kT,  and  ipQ  and  tpQ 
were  calculated  from  Eqs.  (20)  and  (21),  using  parametric 
values  for  the  GaAs-like  and  the  A£As-like  branches  respec- 
tively. 

Ionized  impurity  scattering  was  represented  by  a 
(22) 

relaxation  time 

xII(e)  « 0.41417  K2  T3/2(m*/m) 1/2Ng1  e3/2  (23) 

where 

Ns  = (N+  + N")  g (24) 

and 

1.297xl014  K (m*/m)T2 

I - *n  ^ 1.  (25) 


this  approximate  approach  being  reasonably  valid  so  long 
as  g > 1 or  2, 

Finally,  Eq . (2)  was  used  to  represent  the  alloy-like 
scattering . 

The  energy  dependent  inverse  relaxation  times  of 
Eqs.  (2),  (18),  (19),  (22)  and  (23)  were  summed  to  give  the 
net  inverse  relaxation  time  l/x(e)^22\  and  the  drift 
mobility  was  calculated  numerically  as 


H , * x ( e)  * — - — t -nr 

Md  m*  v 7 m*  .«>  3/2 


/“e3^2x(e) e"£de 


V 


e~ede 


The  Hall  factor  r^  was  calculated  as 


rH  - <T2/-0)/CT»  + i) 


(26) 


(27) 


where 

T1  - T(e)/(1  + <A(e)2 


(28) 


14) 


and 

T2  * woT(e)2/(l+to^T(e)2)  (29) 

and 

_ TJ 

Wq  = ^ is  the  cyclotron  frequency. 

In  Eqs.  (28)  and  (29),  the  bar  represents  the  energy 
weighted  average  over  a Boltzmann  distribution,  as  expanded 
in  Eq.  (26)  for  x (e)  . 

2 

Thus,  given  x,  T,  B,  Ng  and  4E  , the  Hall  mobility 
could  be  calculated  as  r^u^  and  compared  with  the  experimental 
value . 

D.  Data  Fitting  Techniques 

The  computer  program  was  first  tested  on  the  Hall  data 

of  the  GaAs  sample  to  see  if  the  values  of  NA  generated 

were  temperature  independent.  NA  was  found  to  increase 

monatonically  with  temperature  from  4.5  x 10^  cm  ^ at  25°K 

to  8.3  x 10^  cm~^  at  110°K,  indicating  that  the  extent  of 

lattice  and/or  alloy-like  scattering  was  being  underrepresented. 

In  our  attempts  to  make  NA  constant,  we  studied  the  effect 

of  varying  several  parameters.  An  excellent  data-fit  below 

55°K  could  be  obtained  by  including  a small  amount  of  space 

2 2 

charge  scattering  (by  taking  x = .0001  and  Eg  * 9 . 1 eV  in 
Eq.  (2),  equivalent  to  NO  = 377  cm"1),  or  comparably  by 

b 

using  E * 8.3  eV  instead  of  7.0  eV.  To  fit  the  data  in  the 
c 

55°K  - 110°K  range,  we  had  to  reduce  the  polar  optical 
relaxation  time  tpQ(e)  by  a factor  ZZZ  * 0.65.  Although 
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this  adjustment  looked  rather  arbitrary,  it  seemed  to  be 
the  only  one  that  reduced  NA  sufficiently  at  higher 
temperatures  to  make  it  (NA)  temperature  independent,  while 
not  affecting  it  (NA)  for  T < 55°K. 

To  fit  the  Ga,  _ Al  As  data,  it  seemed  appropriate  to 
keep  5C  constant  at  7.0  eV  and  adjust  EBSQ,  since  EBSQ  was 
a variable  parameter  anyway.  Also,  alloy- like  scattering 
was  expected  to  dominate  over  both  D.P.  and  P.E.  scattering 
in  the  temperature  range  studied,  for  x > 0.05  (see  Discussion). 
Hence  the  effect  on  NA  of  adjusting  5 or  h^  would  be 
relatively  small  compared  to  the  effect  of  adjusting  EBSQ. 

Thus  EBSQ  was  adjusted  to  fit  the  data  for  T < 55°K.  Then 

the  empirically  established  procedure  of  adjusting  the 

factor  ZZZ  was  utilized  to  equalize  NA  at  the  higher  temperatures. 

IV.  RESULTS  AND  DISCUSSION 
A.  Results 

The  values  of  EBSQ  and  ZZZ  that  gave  the  best  data  fit 

or  the  near-best  fit  for  the  various  samples  are  listed  in 

Table  B,  along  with  the  values  of  NA  generated.  The 

quality  of  data  fits,  rated  excellent,  good,  fair  or  poor, 

2 

is  also  listed.  Some  values  of  Eg  differ  from  the  corresponding 

EBSQ  by  .01  because  equally  good  "best"  fits  were  obtained 

for  two  adjacent  values  of  EBSQ  (differing  by  .02),  one  of 

2 

which  is  listed.  The  range  of  Eg  indicates  the  range  of 
EBSQ  values  for  which  comparable  constancy  was  obtained  for 
NA.  For  the  GaAs  sample  F-25,  a small  but  finite  value 
of  x (=  .0001)  was  used  to  allow  the  inclusion  of  space 


$• 
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charge  scattering.  The  data  for  sample  F-ll  was  not 
analysed. 

2 2 2 

The  values  of  Eg  obtained  ranged  from  .13  eV  to  .26  eV  , 

corresponding  to  Eg  = .36  eV  to  .51  eV,  the  average  being 
Eg  = .44  eV.  Fig.  3 shows  the  equivalent  values  of  NsQ 
plotted  against  x,  if  the  alloy-like  scattering  is  interpreted 
as  consisting  of  only  space  charge  scattering.  While  the 
data  is  insufficient,  the  extent  of  alloy-like  scattering 
seems  to  be  lower  in  type  (odd  numbered)  than  in  type- 

A (even  numbered)  samples.  The  possible  significance  of  this 
is  discussed  in  Section  IV. C. 

In  addition,  the  compensation  ratio  was  measured  to 
be  about  2 for  most  of  the  samples,  as  shown  in  Table  B. 

B . Sources  of  Error 

The  results  of  this  experimental  analysis  are  subject 
to  uncertainties  and  errors,  both  random  and  systematic, 
from  a variety  of  sources,  namely  data  measurement,  parameter 
estimation,  approximations  made  in  the  mathematial  formulation, 
and  the  data  fitting  technique  for  which  no  a priori  justifi- 
cation could  be  provided.  For  the  apparatus  used,  we  would 
estimate  the  uncertainties  in  B,  y^,  Q^,  T,  and  d to  be 
about  +37.,  +47.,  +107,,  +57,  and  + 0.75y  respectively,  and  the 
uncertainty  in  x for  most  samples  to  be  + .01. 

The  'formulation'  errors  arise  from  the  various  approxi- 
mations made.  For  instance,  the  assumption  1/x  =*  I1/t.  is 

i 3- 

not  strictly  accurate^^.  The  non- degenerate  (Boltzman's) 
approximation  is  largely  valid  for  the  high  purity  samples 
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but  is  in  error  by  a few  percent,  specially  at  the  lower 

temperatures.  The  energy  averaging  of  the  screening  factor 

g in  Eqs.  (24)  and  (25)  instead  of  using  the  energy  dependent 
(22) 

gv  ' can  cause  significant  error,  specially  at  higher  doping 
and  low  temperatures.  Again,  for  n below  1 x 10^,  this  error 
is  small  at  25°K.  Finally  the  dependence  of  space 

charge  scattering  assumed  as  mentioned  earlier,  might  not 
be  correct. 

The  seemingly  arbitrary  need  to  reduce  Tpo(e)  of  Eq. 

(22)  by  a factor  ZZZ  (=.6)  to  fit  the  data  could  be  a result 

of  inaccuracies  in  Eq.  (21)  for  ypQ.  The  solution  of 

Boltzmann's  equation  for  polar  optical  scattering  is  complex, 

(27  30  31) 

and  has  led  to  a variety  of  different  expressions.  ’ ’ ' 
The  term  x(z)  in  particular,  which  we  have  approximated  by 
3/ifz/8  , can  be  in  significant  error  by  a temperature  de- 
pendent factor  which  is  compensated  in  our  analysis  by  ZZZ. 
Compounding  this  formulation  error  is  the  high  sensitivity 
of  the  calculated  value  of  ypQ  to  some  of  its  parameters. 

For  instance,  a 17.  error  in  estimating  either  Kg  or  Kq  would 
lead  to  a 77.  error  in  ypQ,  as  would  a 17.  error  in  z at 
60°K.  Finally,  using  a relaxation  time  approximation  for 
ypQ  itself  could  have  caused  some  error  which  required 
compensation. 

Finally,  to  illustrate  the  effects  of  errors  in 
estimating  the  various  physical  parameters,  we  select 
sample  F-17  (x  - .11).  Fig.  4 shows  the  calculated  values 
of  yp£,  yDp,  ypo,  yAL  and  yIX  plotted  against  T.  The 


mobilities  correspond  to  the  estimated  or  measured  values  of 
the  relevant  parameter  for  sample  F-17,  except  for  UpQ»  which 
is  shown  for  ZZZ  ■ 1.  For  calculating  Upj.  we  assumed 

nd  - 3na- 

It  is  seen  that  between  25°K  and  100°K,  ionized  impurity 
scattering  and  alloy-like  scattering  dominate  over  the  others, 
which  implies  that  the  effect  of  varying  upE,  HDp>  and  Up0 
on  Vtota].  is  small,  especially  in  the  lower  half  of  the 
temperature  range.  For  example  a 107=,  error  in  estimating 
h^  would  cause  only  a 1%  error  in  Eg.  However,  for  samples 
with  lower  x,  y^is  higher  and  the  value  of  Eg  is  more 
susceptible  to  errors  in  the  estimation  of  lattice  scattering. 
For  temperatures  below  ~60°K,  ionized  impurity  scattering 
dominates,  so  that  the  values  of  NA  generated  are  more  accurate. 
At  higher  temperatures,  Ng  + is  estimated  from  a smaller 
component  of  the  total  scattering,  and  is  therefore  subject 
to  larger  systematic  errors,  which  might  also  contribute  to 
the  ZZZ  factor. 

C . Discussion 

Both  Kaneko  et  al.^^  and  Stringfellow^^  report  a 
good  theoretical  fit  to  their  y vs  T data  by  assuming  the 
space  charge  scattering  parameter  NgQ  to  vary  as 

N Q - 5xl03  + 6 . 3x105x  . (30) 

s 

By  substituting  Eq.  (30)  in  Eq.  (4)  and  comparing  with  Eq . (2) 
at  x * .1,  .2  and  .3,  it  is  seen  that  space  charge  scattering 
described  by  Eq.  (30)  is  equivalent  to  alloy-like  scattering 
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2 2 2 
with  Eg  = 1.64  eV  . In  constrast,  we  obtain  values  of  Eg 

2 

from  .13  to  . 26  eV  , or  equivalently 


NgQ  * 8.6x10  x 


(31) 


I 


A 


r 


The  samples  from  references  (9)  and  (10)  were  doped 
17  -3 

in  the  range  of  10  cm  , about  two  orders  of  magnitude 

higher  than  our  samples.  Given  this  limited  information, 

some  correlation  is  seen  to  exist  between  the  extent  of 

space  charge  scattering  and  the  sample  purity.  Even  among 

2 

our  samples,  we  measure  lower  values  of  Eg  for  the  type  B 

samples  compared  with  the  less  pure  type  A samples , although 

layer  F-19  (x  * .154,  NsQ  = 16600)  violates  this  rule. 

We  find  it  difficult  to  draw  conclusions  about  the 

extent  of  alloy  scattering  in  our  samples , except  that  AE 

is  definitely  lower  than  AEg  - .7  eV.  The  fact  that  the 

2 

measured  values  of  Eg  differ  by  as  much  as  2:1  suggests  the 

2 

presence  of  some  space  charge  contribution  to  Eg,  though 
we  cannot  rule  out  the  possibility  that  this  variation  re- 
sults mainly  from  measurement  and  analysis  errors.  If 
Ferry' s^  estimation  that  AE  = .12  eV  is  indeed  correct, 

then  the  alloy  scattering  contribution  to  the  measured  value 
2 

of  Eg  will  be  negligible.  On  the  other  hand,  if  the  value 

(3} 

of  Ae  predicted  by  Harrison  and  Hauser'sv  1 model,  (55eV) , 

is  reasonably  true,  then  alloy  scattering  will  contribute 

2 

a major  part  of  Eg. 


D.  Summary 

Hall  mobilities  of  high  purity  LPE  Ga, AH  As  layers 

“ X X 

were  measured  as  a function  of  temperature  in  the  range 
( 0 _<  x < .2).  The  various  physical  parameters  that  are 
relevant  to  electron  scattering  were  estimated  for  Ga,  AH  As , 

1 " X X 

using  the  best  extrapolation  techniques  available  to  us. 

A computer  program  was  used  to  generate  Hall  mobilities  as 
a function  of  the  independent  experimental  and  sanple  parameters, 
namely,  x,  T,  B,  and  Ng , and  to  adjust  Ng  for  test  values  of 
Eg  to  obtain  agreement  with  experiment.  From  Ng  and  the 
carrier  concentration,  the  acceptor  density  NA  was  calculated, 
and  its  constancy  with  temperature  was  used  as  the  criterion 
for  selecting  Eg  from  its  test  values.  The  alloy-like 
scattering  parameter  Eg  was  measured  to  lie  between  . 36  eV 
and  .51  eV,  the  average  value  being  .44  eV.  The  extent  of 
alloy-like  scattering  in  our  samples  was  thus  an  order  of 
magnitude  lower  than  in  samples  of  references  (9)  and  (10). 

No  definite  conclusions  could  be  reached  about  the  value  of 
AE,  the  depth  of  the  random  potential  well  in  Ga,  AH  As , 

.L  “ X X 

except  for  placing  an  upper  limit  of  ~ . 5 eV  on  its  value. 
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Table  A 

Thickness  and  Composition  Data  for  the 
Ga-  AH  As  (0<x<.18)  Studied 


Sample 

# 

Thickness 

(microns) 

Photolum-  Q 

inescence  (5  K) 
peak  energy 
(eV) 

x,  from 

p . 1 . 

% 

x from 
thermo- 
dynamics 

% 

F-25 

14.0 

0 

F-10 

n.i 

1.545 

2.0 

2.1 

F-ll 

ii.i 

1.541 

1.75 

1.75 

F-12 

12.0 

1.581 

4.6 

5.1 

F-13 

13.9 

1.577 

4.2 

4.3 

F-14 

11.4 

1.632 

8.2 

8.8 

F-15 

13.0 

1.625 

7.8 

7.7 

F- 16 

11.1 

1.682 

11.8 

12.6 

F-17 

10.3 

1.670 

11.0 

11.3 

F-18 

9.3 

1.744 

17.7 

17.1 

F-19 

9.5 

1.726 

15.4 

15.6 
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FIGURE  CAPTIONS 

1.  Hall  mobility  vs  temperature  data  for  Ga..  Al  As  Samples. 

1“  X X 

2.  Hall  carrier  density  per  unit  area  (QH)  3 temperature 
data  for  Ga1  Al  As  samples.  The  discontinuities  were 
caused  by  a 'set-up'  error. 

3.  Space  charge  parameter  NsQ  (assuming  all  alloy-like 

scattering  to  be  due  to  space  charge  scattering)  plotted 

against  x.  The  vertical  lines  are  error  bars  estimated 

2 

from  the  uncertainty  in  measuring  Eg. 

4.  The  calculated  temperature  variation  of  the  various 

mobility  components  for  representative  sample  F-17 

(x  = .11).  Values  of  the  physical  parameters  used  in 

7 -1 

the  calculations  are:  5c  = 7.0  eV,  h^4  = 1.4x10  V • cm  , 

E2  = .15  eV2,  N_  + N.  = lxlO15  cm-3  and  = 3N..  For 
B 'DA  DA 


the  UpQ  curve,  ZZZ  was  assumed  to  be  unity.  The  curve 
marked  y represents  the  experimental  Hall  mobility. 
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